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ABSTRACT 

Context. Despite many studies of star formation in spiral galaxies, a complete and coherent understanding of the physical processes 
that regulate the birth of stars has not yet been achieved, nor has unanimous consensus been reached, despite the many attempts, on 
the effects of the environment on the star formation in galaxy members of rich clusters. 

Aims. We focus on the local and global Schmidt law and we investigate how cluster galaxies have their star formation activity 
perturbed. 

Methods. We collect multifrequency imaging for a sample of spiral galaxies, members of the Virgo cluster and of the local field; we 
compute the surface density profiles for the young and for the bulk of the stellar components, for the molecular and for the atomic 
gas. 

Results. Our analysis shows that the bulk of the star formation correlates with the molecular gas, but the atomic gas is important 
or even crucial in supporting the star formation activity in the outer part of the disks. Moreover, we show that cluster members that 
suffer from a moderate HI removal have their molecular component and their SFR quenched, while highly perturbed galaxies show 
an additional truncation in their star forming disks. 

Conclusions. Our results are consistent with a model in which the atomic hydrogen is the fundamental fuel for the star formation, 
either directly or indirectly through the molecular phase; therefore galaxies whose HI reservoirs have been depleted suffer from 
starvation or even from truncation of their star formation activity. 
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1. Introduction 

A satisfactory understanding of the physical processes governing 
the formation of giant molecular clouds (GMC) from primeval 
atomic hydrogen (HI) and their instability leading to the forma- 
tion of stars is still far from achieved. This is not surprising due 
to the complex mix of hydrodynamical and gravitational physics 
involved in the process of star formation. Beside the theoretical 
difficulties, what surprisingly limits our present understanding of 
this issue is the lack of a robust phenomenology of star formation 
in nearby galaxies. Despite the many recipes that have been pro- 
posecQ, the simples t and most com monly used parametrization 
is the Schmidt law (Schmidt 1959) that establishes a correlation 
between the star formation rate density p s f r and the gas density 

Pgas' 

Psfr = ap a gas ; (1) 

derived for our Galaxy, the v alidity of this la w has been tested 
also for external galaxies (e.g. lKennicuttl 19831) . for which a sim- 
ilar relation, usually written in terms of surface densities, holds: 

^sfr=AJT gas . (2) 

Many studies of the global Schmidt law, i.e. the relationship be- 
tween the disk-averaged star formation rate (SFR) and the total 



1 See for i nstanc e the dynamically modified Schmidt law in 
I Wong & Blitz! d2002l) 



gas content, have been carried out on l arge samples of g alax- 
ies. A milestone paper on this issue is bv lKennicuttl (fl998). who 
studied the connection between the integrated SFR surface den- 
sity and the integrated gas surface density in a sample of normal 
and starburst galaxies. Using integrated quantities, he found that 
the same Schmidt law with an index n — 1 .4 holds for both nor- 
mal and starburst galaxies, covering a range of 5 dex in the gas 
density and over 6 dex in the SFR; he also found that this re- 
lation is mainly driven by the correlation between the SFR and 
the HI, while there is a poorer correlation between the SFR and 
the H2, significant only when massive (L^ > 10 10 L Q ) objects 
are considered. Owing to the fact that star formation takes place 
in molecular regions on parsec scales, this result was somewhat 
unexpected, but Kennicutt argued that the poor knowledge of 
the CO-to-H2 conversion factor might be responsible for the 
scatter in the H 2/SFR correlation. Using seven galaxies fro m the 
BIMA survey dHelfer et alj 120031) . IWong & Blitzl d2002l) stud- 
ied the applicability of the Schmidt law on local scales, i.e. 
through a comparison between the SFR and the gas surface den- 
sity profiles. They found that a Schmidt law between the SFR 
and the gas content holds also on local scales and, even if the 
molecular gas content alone does not always control the SFR, 
they concluded that the correlation found for the total gas is 
entirely driven by the molecu lar component. Consistent results 
were found by iBoissier et al.l (120031) . who used low resolution 
CO data but included an X factor that varies as a function of 
the metallicity. Another open issue is to asses the influence of 
the environment on the star formation. Galaxies in clusters suf- 



2 



Fumagalli & Gavazzi: Gas content and star formation rate in spiral galaxies. 



fer from a variety of environmental perturbations, mainly due 
to gravitational interactions between galaxies themselves or be- 
tw een a galaxy and the clu ster potential well (see the review 
by Boselli & Gavazzi 2006) or to various hydrodynamical inter- 
actions between the ISM and the intergalactic medium (IGM). 
Tidal interactions between gala xies or betwee n individual galax- 
ies an d the cluster potential dMerrittl Il984t iBvrd & Valtonenl 
1990) produce both the removal of the outer and looser com- 
ponen ts and the gas inf all towards the galaxy center; harass- 
ment (lMooreetal.| [l996) can induce sinking of gas towards the 
galaxy ce nter and can shape th e stellar profiles; the ram pr essure 
stripp ing dGunn & Gottlll972l) or the viscous stripping (Nulsen 
1982) can effectively remove the gas components, mainly from 
the outer part of the disks. Furthermore, a combination of gravi- 
tational and hy drodynamic effects can cause ga laxy starvation or 
strangulation dLarson et al.l fl980; Bekk i et alll2002l) . According 
to starvation, the gas that feeds the star formation in the lo- 
cal universe probably comes from the infall of an extended gas 
reservoir, therefore the effect of removing the outer galaxy halo 
would be that of preventing further gas infall; as a consequence, 
the star formation exhausts the available gas, quenching fur- 
ther activity. The literature about the classification and the ef- 
fects produced by the envir onment on the evolutio n of clus- 
ter galaxies is very rich (e.g. iBoselli & Gavazzil 2006), but, de- 
spite many studies on nearby and distant clusters, it is not clear 
what the dominant processes are, nor how the star formation 
activity is affected (Koopma nn & Kennevl 12004. and reference 
therein). Unanimous consent seems only to hold on the statement 
that environmental processes are effective at remo ving the outer 
ISM and therefore affecting the HI component dCavatte et alj 
Il994t iGiovaneili & Havneslll985l) . leaving the molecular com- 
ponent, well-bounded inside the galaxy potential w ell, undis- 
turbed dKennev & Young|ll989t IBoselli et alJl2002bl) . In order 
to study the environmental effects on star formation, we col- 
lect both objects that are members of the Virgo cluster and lo- 
cal field galaxies. Since the different components involved in the 
star formation are observable along a broad stretch of the elec- 
tromagnetic spectrum, we adopt a multifrequency analysis. To 
quantify the atomic hydrogen we collect observations at 21 cm; 
the molecular hydrogen content is estimated indirectly via CO 
observations at 2 mm; we quantify the stellar mass with obser- 
vations at the near-infrared (NIR) or visible bands and we study 
the presence of new stars indirectly, observing the hydrogen re- 
combination line (Ha) at 6563 A. The sample is presented in 
Sec. [2] together with the description of the data reduction proce- 
dures; the analysis is given in Sec. [3] while we discuss our results 
and conclusions in Sec.[4]and[5] 



2. Data 

2.1. The sample 

The sample of galaxies chosen for the present analysis is se- 
lected primarily according to the availability of high sensitivity, 
extensive CO mapping, s uch as in th e Nobe yama CO atlas of 
nearby spiral galaxies by iKuno et alJ (l2007h . Secondly we re- 
quire that reliable Ho- imaging exists in the literature. 28 giant 
spiral galaxies match these criteria. For all of them we were able 
to find additional imaging material in the red continuum (i or H 
band) and for 25/28 we found sensitive HI mapping. Table[T]lists 
some observational parameters and references for the selected 
galaxies. The sample is well balanced between isolated objects 
and members of the Virgo cluster, as stressed by the local den- 
sity parameter (see Col. 1 1). One characteristic of our sample is 



that it is composed of giant spirals with a typical H-band lumi- 
nosity of 10 10 - 10 1 1 L . No BCDs nor Irr/dlrr are included. The 
fact that galaxies in the Nobeyama atl as have not been selected 
according to a CO emission criterion dKuno et al.ll2007l Sec. 2) 
protects against a possible bias towards CO rich galaxies. One 
could question why we preferred the Nobeyama atlas to higher 
resolution interfero metric CO material, such as the BIMA Survey 
of Nearby Galaxies dHelfer et al. 2003) that provides images for 
44 nearby spiral galaxies. This is mainly to ensure that our flux 
measurements are unaffected by severe missing flux problems on 
scales larger than 20-30" at 1 15 GHz. Even though spectra from 
the central galaxy regions have been used to recover the missing 
flux or single dish observations provide the zero-spacing data, 
the algorithms used to recover th e extended flux are somewhat 
unreliable (see lHelfer et al .1120031 App. A). 

2.2. Ha data 

We compute the star formation rate fro m the intensity of the Ha 
emission line, following the method o f lKennicuttldl983L [i998). 
Ha images for the V CC galaxies are retrieved from GOLDMine 
dGavazzi et alj 120031) . while for the remaining galaxies we re- 
trieved images fro m NEE0- Images f r om G OLDMi ne are cal- 
ibrated following [Boselli & Gavazzi (2002]) and IBoselli et al.l 
(2002a), i.e. using the photometric zero-point included in 
the image headers. For the images collected from NED in- 
stead, we derive the photometric calibration using published 
fluxes dYoungetal.lll996t iRomanishinll 19901: iKennicutt & Kent! 
|1983|P 1 Before converting the Ha flux into a star formation rate, 
we correct for the [Nil] emission; for almost all the galaxies we 
apply a correction derived from drift-scan spectra that accounts 
for the mean [Nil] emission from the entire disk. Where individ- 
ual correction factors are not available, w e apply the standard 
value [Nil] /Ha = 0.53 dKennicuttl 1 1 9921) . Since we consider 
[NII]/Ha factors averaged over the entire disk, the correction un- 
derestimates the [Nil] emission from the bulges of galaxies that 
host active galactic nuclei (AGN), where nuclear [Nil] can be 
even stronger than Ha. This is particularly severe in our sample 
composed of giant spiral s, since almost all m assive galaxies are 
known to host an AGN dDecarli et al.ll2007l) . In addition, even 
if the nuclear emission was corrected by the suitable factor, it 
would not be possible to disentangle the Ha flux due to star for- 
mation from the one due to the AGN ionization. Therefore, by 
inspecting both the Ha imaging and the available spectra, we 
identify a nuclear region for each galaxy, and in the following 
analysis we do not include the emission from these regions. As 
for the dust absorption, since the Ha/H/3 or the Ha/Br y ratios 
are unknown for the majority of our galaxies, we prefer to ap- 
ply a standard dust extinction A = 1.1 mag dKennicutti ri983) 
to other more sophisticated, yet unreliable corrections. From the 
corrected Ha lumi nosity we comp ute the star formation rate us- 
ing the conversion dKennicuttl 1998b 



SFR [Mayr 1 ] = 7.94 X 10~ 42 L(Ha) [erg/s] 



(3) 



based on a Salpeter IMF (a = 2.35) over 0.1 - 1OOM . 
Furthermore we assume that a fraction /=0.57 of Lyman contin- 
uum photons ionizes the hydrogen ato ms inside the star for ming 
regions, where dust is mixed with gas (Hirash ita et al ]|2003l) . On 
the contrary, we neglect the correction for the escape fraction 
of continuum Lyman photons, since it has been shown that it is 



2 See table[T]for individual references. 

3 We correct fluxes from Kennicu tt & Kentl £l983) with the factor 
1.16, as suggested by the authors 
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Table 1. Observational parameters. Column 1,2: NGC, VCC or UGC designations. Column 3,4: coordinates (J2000) from 
GOLDMine and HyperLeda. Column 5: recessional velocities in km/s from GOLDM ine and HyperLeda. Column 6: distances 
in Mpc, a ssign ed following the subcluster membership criteria of Gav azzi et alj d 19991) or from the Nearby Galaxies Catalogue 
by iTullyl (Il994l) . Column 7: galaxy inclinations in deg from GOLDMine and HyperLeda. Column 8: morphological types from 
GOLDMine and HyperLeda. Column 9: H-band luminosity in L from GOLDMine or computed from H magnitude from SIMBAD. 
Column 10: radius in arcsec at the B-band 25th isophote from HyperLeda. Column 11: density i n Mpc' 3 f or galaxies brighter than 
-16 mag in the vicin ity of the galaxy consider ed, de termined on a 3D-grid a t 0.5 Mpc spaci ng (Tullv 19941) . Column 12: referenc e 
to Ha imaging: a = lBoselli & Gavazzil d2002l); b = iKoopmann et al l (120011); c = iBoselli et all d2002al): d = iMarcum et ail d2001l) ; 
e - iKnapen et alj d2004l): f = iMeurer et al.l d20()6l) : g = IKennicutt et all d2003l): h = iJames et all (12004b . Column 13: reference to 
metall icity: i = ISkillman et al.1 d!996l) : i = GOLDMine: k = IKennicutt et al l d2003l): 1 - Ivan Zee et alj (fl998h: m = IZaritsky et all 
d!994l) . Column 14: reference to HI mapping: n = THINGS; o = VIVA; p = ICavatte etail d 19901): a = IVerheiien & Sancisil d2001l) . 
Column 15: reference to NIR/Visual imaging: r = SDSS; s = GOLDMine; t = lKuchinski et all d2000h ; u = lMulchaev et all d 19971) . 



less than 6% in starburst galax ies and even negligible for normal 
galaxies (Hec kman et alj200lb . We estimate that the typical flux 
unce rtainty is 0.05 mag for V CC galaxies; 0.1 mag for fluxes 
from IKennicutt & Kent! (TT9831), 5% for those from I Young et al.l 
d 19961) and up to 20% for lRomanishir] d!990l) . These errors pro- 
duce a final uncertainty on the star formation rate which ranges 
from 0.04 dex to 0.09 dex, with a mean value of 0.05 dex. An 
additional uncertainty of 0.2 dex affects all data due to the un- 
certainty in the dust correction. 

2.3. H 2 data 

In the absence of a direct tracer of the H2 molecule, the molec- 
ular hydrogen abundance can be inferred from CO observations. 
We collect flux calibrated maps of the I2 CO(J : 1 -0) transition at 
115 GHz from the Nobeyama CO atlas o f nearby spira l galax- 
ies taken with the 45-m radio telescope by iKuno et ail d2007l) . 
Despite their spatial resolution of 15 arcsec (at 115 GHz), these 
observations have been preferred to higher resolution interfero- 
metric data because of their higher sensitivity to the large-scale 



emission from galaxy disks. In fact the typical rms noise of 40- 
100 mK combined with a velocity resolution of 5 km/s, gives a 
sensitivity of N(H2) — 5 - 10 x 10 19 cm _1 . The most challenging 
task related with the CO observations, both from a theoretical 
and an observational point of view, is the determination of the 
CO-to-H2 conversion factor (X in cirT 2 /(K km/s)). Many at- 
tempts to determine a reliable CO-to-H2 conversion factor can 
be found in the literature; nevertheless, there is not yet unan- 
imous consent on the dependence of X on parameters such as 
density, temperature or even th e geometry of the clumps that 
constitute the molecular regions (Wall 2007). Due to these diffi- 
culties, a constant conversion factor X, as derived for our Galaxy, 
has been often assumed, probably underestimating the H2 in low 
luminosity galaxies. In addition, it is expected that the CO-to- 
H2 conversion strongly depends on dust and metallicity that reg- 
ulate the H2 formation and photodissociation or the CO heating. 
Therefore, we avoid a single conversion factor for our sample, 
but we choose individual X as a fu nction of the metallic ity, fol- 
lowing the empirical calibration bv lBoselli et al.l d2002bl) . made 



4 



Fumagalli & Gavazzi: Gas content and star formation rate in spiral galaxies. 





50 100 150 300 

Radius (arcsec) 



Fig. 1. Example of multiband images and profiles: NGC 4051. Top-left: i band; top-right: Ha; bottom-left: CO; bottom-right: 
surface brightness profiles. Y, s f r (red) in M Q pc~ 2 Gyr -1 ; I. HI (magenta), (blue) and l. s , ar (green) in M pc~ 2 . The arrow indicates 
r 25 as in Col. 10, Tab. Q] 



on fourteen independent measurements of CO and molecular hy- 
drogen. 

logX = -1.01 + 0.14 x (12 + log(0/H)) + 29.28 + 0.20 . (4) 

With this equation we compute the conversion factor for each 
galaxies (See Col. 13 in Table [T] for references to metallicity 
data). For seven galaxies without published 12 +log(Q/H) we de- 
rive the metallicity using (Boselli et al. 2002b|) 

12 + log(0/H) = 0.37 + 0.04 x log Lh [L ] + 4.98 + 0.15 , (5) 

calibrated on a sample of nearby and Virgo cluster galaxies. 
The mean discrepancy between the metallicity directly measured 
from spectra and derived using © is 0.2 dex in our sample, not 
much greater than the uncertainties on the measurements them- 
selves and significantly lower than the statistical uncertainty on 
the X conversion factor derived using (0J, that is on average 1 .3 
dex. One should keep in mind that other unknown dependen- 
cies, beside the one on metallicity, would increase the uncer- 
tainty on X. Since spiral galaxies exhibit a radial metallicity gra- 
dient, when available in the literature, we account also for them 



to compute X as a function of the galaxy radius. We find that 
X spans 0.4 dex on average in individual galaxies and in some 
cases up to 1 dex. With the X conversion factor, we compute the 
molecular hydrogen column density using: 

N(H 2 ) [M cm" 2 ] = 2 x 1 .36 x m p x X x I co , (6) 

where 1.36 is the helium correction factor, m p is the proton mass 
and I co is the CO flux in K km/s. 

2.4. HI data 

To study the distribution of the HI mass, we collect images at 
21 cm from various sources^. For seven Virgo galaxies, we ob- 
tained moment maps from VIVA, The VIA Ima ging of Virgo 
in Ato mic gas (Chung et al. 2008, in preparation; Chun getail 
2007); VIVA images have a spatial resolution of 15 - 17 arc- 
sec and they reach a typical sensitivity of 1.5 x 10 19 crrT 2 . For 

4 See Table[T]for individual references. 
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nine field galaxies, instead, we obt ain HI images from THINGS, 
The HI Nearby Galaxy Survey dWalter et alj [2008h ; THINGS 
images have a spatial resolution of 11 arcsec and a sensitiv- 
ity of 1.2 x 10 19 cm 2 . For the remain ing galaxies, we coll ect 
HI d ata from the literature, eith er from ICavatte et alj d 19901) or 
from Yerheijen & Sancisil d200ll) . These maps have a lower spa- 
tial resolution, spanning from 40 to 70 arcsec. The sensitivity 
ranges from 1 x 10 19 cirT 2 to 1 x 10 20 crrT 2 . For three galax- 
ies (NGC 1068, NGC 4212 and NGC 6951) we cannot find HI 
maps, but only a measurement of their integrated HI mass. From 
the flux calibrated images, we compute the neutral hydrogen sur- 
face density using 



N(HI) [M cm- 2 ] = 1.36 x m p x l m . 



(7) 



where 1.36 is the correction factor for the helium and Ihi is the 
column density in cirT 2 . The total HI mass in solar units is com- 
puted using: 



logM ffl [Mq] = log(2.36 x 10 5 x SintP x D 2 ) , 



(8) 



where the total HI flux SintP (Jy km /s) is from ALFALFA The 
Arecibo Legacy Fast ALFA Survey dGiovanelli et al.1 120051) . or 
from the HyperLeda database and D is the galaxy distance in 
Mpc. 

2.5. Optical data 

The stellar mass distribution is obtained from imaging in 
the NIR H band from GOLDMine or the i band from the 
Sloan Digital Sky Survey (SDSS-DR6. lAdelman-McCarfhv et all 
2008). Exceptions are NGC 5236 and NGC 6951 (not imaged 
in the SDSS) for which we collect respectively I band and K 
band images from NED. In Table Q] we show individual refer- 
ences. Due to their large angular dimension, our galaxies are 
often cut in two or more SDSS fields, thus we retrieve im- 
ages using Montage, a toolkit for assembling FITS images over 
multiple frames into a mosaic. The photometric calibration for 
galaxies from GOLDMine relies on the zero-points provided i n 
the image headers, obtained according to lGavazzi et al.l (fl996b); 
GOLDMine images for NG C 4192 and NGC 456 9 are from The 
Two Micron All Sky Survey (Skru tskie et al.l l2006), thus for these 
galaxies we adopt the 2MASS standard calibration procedure. 
For all images, the typical uncertainty is <0.05 mag. SDSS im- 
ages are calibrated (within 1% accuracy) using the zero-points 
derived from a set of parameters (calibration constant, air mass 
and extinction coefficient) provided in the tables associated with 
the images. To derive the stellar mass, we apply to each galaxy 
a mass-to-light ratio , computed throu gh an extrapolation along 
the IMF according to lBell et all (120031) : 



logCM/L),, [Mo/Lq] = a A + ■ Color , 



(9) 



with ax and b\ chosen accordingly to the band used. We use B- 
V Color from HyperLeda in combination with NIR data and g—i 
for i band data. The mean dispersion on M/L is 0. 1 dex. 

2.6. The surface brightness profiles 

For each galaxy we derive a calibrated surface brightness pro- 
file at four frequencies. All images are firstly astrometricalry 
calibrated using GAIA, then aligned and resampled to 3 arc- 
sec/pixel. Stars near the galaxies are masked using the IRAF 
task imedit and the local sky background is estimated and 



subtracted using marksky and skyfit from Galphofl Using a 
Gaussian kernel we smooth all images to a spatial resolution of 
15 arcsec, in common with the CO observations and most of the 
HI images. To derive the surface brightness profiles we use the 
IRAF task ellipse that fits elliptical isophotes to the galaxy 
images. Starting from a set of initial parameters given manu- 
ally for the NIR or / band imaging, the fit maintains as free pa- 
rameters the ellipse center, the ellipticity and the position angle; 
the ellipse semi-major axis is incremented by a fixed fraction 
of its value at each step of the fitting procedure and the rou- 
tine halts when the surface brightness found in a given annulus 
equals twice the sky rms. This task provides the statistical er- 
ror of the surface brightness. From the same set of ellipses we 
derive the surface brightness profiles f or the remaining im ages, 
with the exception o f the H I data from Cav atte et al.l (1 19941) and 
IVerheiien & Sancisil d2001l) for which we interpolate the pub- 
lished profiles on our grid. The derived surface density profiles 
are corrected for projection effects using the galaxy inclination 
on the plane of the sky. For each galaxy we finally obtain the 
HI, H2 and stellar surface density profiles E# 2 and S Jfar 
in M Q pc~ 2 ; the star formation rate surface density profiles 2 s / r 
are in M Q pc~ 2 Gyr -1 . Moreover, we compute the color profiles: 
sfr-sta = logE s/r -logS Jfflr , sfr-H 2 = logE j/r -logS ff2 and 
H2 - sta = log "Eh-, -log£ J/fl ,-. Fig.[T]illustrates an example of the 
imaging material and density profiles for the galaxy NGC 405 1 . 
From top-right (anticlockwise): Ho- net; /-band (SDSS); CO im- 
age; surface brightness profiles of stars (green), HI (magenta), 
SFR (red) and H 2 (blue). 



3. Analysis 

3.1. The Schmidt law 



Following [Kennicuttl dl998l) . IWong & Blitzl d2002l) and 



Boissier et al. (2003]) we test the validity of the Schmidt law on 
both local and global scales. 

3.1 .1 . The local Schmidt law 

Figure [2] (left) shows the correlation between the SFR and the 
H2 surface brightness in the individual galaxies belonging to the 
studied sample; in doing this comparison, we exclude the nu- 
clear regions as previously discussed. Since in each galaxy the 
profiles do not have exactly the same radial extent, we prefer to 
truncate the most extended profile at the radius of the shortest 
one, avoiding any extrapolation. Typically, the SFR and the H2 
emission have a comparable spatial extent, in which on average 
there is a fair correlation between the SFR and the molecular hy- 
drogen surface density. In fact 80% of our sample shows a cor- 
relation index r > 0.8. Only a few galaxies show a weak correla- 
tion, mainly due to peculiarities such as prominent bars empty of 
star formation (NGC 4303), starburst rings (NGC 4321) or sig- 
nificant gas displacements due to known severe environmental 
perturbations (NGC 4548). For the correlated galaxies, we per- 
form a least square fit; despite the overall validity of the Schmidt 
law, 8/20 objects cannot be fitted using a single power law or 
the obtained rms is greater than 0.15 dex, arbitrarily assumed 
as a threshold. Inspection of the left panel of Figure [2] shows 
that the faintest points in the SFR/H2 relation that belong to the 
outer parts of the disks significantly deviate from the correlation 



5 GALPHOT was developed for IRAF-STSDAS by W. Freudling, J. 
Salzer, and M . P. Ha ynes and adapted by S.Zibetti and L. Cortese. See 



iGavazzietall feOOl). 



6 



Fumagalli & Gavazzi: Gas content and star formation rate in spiral galaxies. 




Fig. 2. The local Schmidt law for the entire sample, separately for the molecular (left panel) and the total gas (right panel). Profiles 
are normalized to the maximum values for each galaxy. The dashed line in the left panel represents the index n = 1, in the right 
panel n - 1 .4. 



driven by the points associated with the inner part of the galaxy 
disks. For the remaining 60% of the correlated galaxies, we ob- 
tain a reliable fit using a single power law; the mean index ob- 
tained is n — . 93 ±0 .11, in good agreement with the finding of 
IWong & B litz (2002) (n = 0.8 + 0.3). Despite this overall agree- 
ment, there are slightly differences in the individual indexes n for 
the five galaxies we have in common, mainly due to the fact that 
we exclude the innermost galaxy regions and that Wong & Blitz 
choose a conversion factor X that does not vary as a function of 
the radius nor from one galaxy to another. 
We checked that there is no evident correlation between the star 
formation rate and the atomic hydrogen, especially in the inner 
part of the galaxy disks where the Ha emission is often strong 
but usually there is a low abundance of HI. To study the role of 
the total gas in the Schmidt law we compute the surface den- 
sity profiles Y,hi+h 2 -> combining the atomic and the molecular 
hydrogen profiles without any extrapolation. In Figure [2] (right 
panel) we plot the Schmidt law for the total gas surface density 
on the same scale used for the molecular hydrogen alone (left 
panel). Figure [2] shows that the addition of the atomic hydro- 
gen makes the correlation between the SFR and the gas better 
than considering the H2 alone; the most striking feature is that 
the faintest points in the local Schmidt law that deviate the most 
when considering the molecular hydrogen alone, now follow a 
single power law. Furthermore, we notice that the dispersion in 
the inner part of the disk is on average reduced, perhaps because 
by adding the atomic component we are taking into account the 
HI produced by the photodissociation in regions that have a high 
star formation activity. Inspecting individual galaxies, we find 
that all but two galaxies with r > 0.8 can be fitted using a single 
power law with an rms < 0.15, i.e. the 60% of the entire sam- 
ple can be described with a single Schmidt law. Using a least 
square fit, we derive a mean index n = 1.39 + 0.10, in good 
agreement with the value obtained by Kennicutt using integrated 
quantities. Despite confirming the correlation between the SFR 
and the H2 that was suggested by Wong & Blitz, the important 
difference between our results and theirs is that the addition of 



HI improves this correlation, significantly reducing the number 
of galaxies for which the Schmidt law fails in the outer part of 
the disks. Therefore we cannot agree with the conclusions of 
Wong & Blitz that the correlation for the total gas is entirely 
driven by the molecular component. A possible cause for this 
disagreement might be their use of CO imaging from the BIMA 
survey. Inspecting the H2 radial profiles for the five galaxies we 
have in common, we detect molecular hydrogen emission that is 
up to two times more extended than in the BIMA observations. 
The analysis of Wong & Blitz seems more limited to the inner 
part of the galaxy disks, for which we agree that the correlation 
is almost entirely driven by the molecular component. Due to 
the superior sensitivity of the Nobeyama mapping to large-scale 
emission, we can push this analysis to the outer parts, where we 
find that the HI is the most important quantity that drives the 
correlation between the star formation rate and the total gas. 

3.1.2. The global Schmidt law 

After having studied the local Schmidt law, we focus on its ap- 
plicability while considering integrated quantities; therefore, we 
compute for each galaxy the total SFR and the gas surface den- 
sity, by integrating along the profiles. In Figure [3] we show the 
correlation between the SFR and the H2 (left panel) and the 
correlation between the SFR and the HI+H2 (right panel), for 
our measurements (filled circ les) and fo r the g alaxies we have 
in common, using data from Kennicutt (1998) (open circles). 
Uncertainties are not drawn, but as previously discussed the SFR 
has a typical error of 0.2 dex while the uncertainty for the H2 
may be as large as 1 dex; the uncertainty for the total gas is 
entirely driven by that of the molecular hydrogen. We see that 
there is a fairly good agreement between our values and the ones 
from Kennicutt, confirming the validity of our calibration proce- 
dures, in spite of the slightly different corrections. Considering 
the left panel, we find a good correlation (r = 0.76) between the 
SFR density and the molecular hydrogen density; therefore, for 
our sample composed of massive galaxies, we can assume that 
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Fig. 3. The Schmidt law for integrated quantities, separately for the mol ecular (left panel) and the total hydrogen (right panel). 
Filled circles are for our measurements, while empty circles are data from lKennicuttl ([1998), for the galaxies we have in common. 
Uncertainties are 0.2 dex for SFR and up to 1 dex for H2 and for the total gas. 



even on a global scale the Schmidt law is a suitable parametriza- 
tion for the relationship between the SFR and the molecular gas. 
Adding the HI (right panel), we can see the same good correla- 
tion (r - 0.77) between the SFR density and the total hydrogen 
density; therefore, considering a sample of massive galaxies, we 
find that the Schmidt law for the total gas does not differ sig- 
nificantly from the one for H2 alone, thus the bulk of the star 
formation that takes place in the inner part of the disks corre- 
lates with the molecular hydrogen surface density, in agreement 
with what is found on local scales. This finding supports the hy- 
pothesis that the poor correlation with the H2 found by Kennicutt 
is driven by the inclusion of low luminosity galaxies for which 
it is difficult to obtain reliable CO observations and for which a 
constant X conversion factor does not hold. 



3.2. The environmental effects 

In the previous section we have shown that the atomic hydro- 
gen seems to play some role in the star formation process, es- 
pecially at the edge of galaxy disks; moreover, the atomic hy- 
drogen, built up during the cosmological nucleosynthesis, is the 
primeval constituent of the molecular hydrogen and therefore 
provides the essential fuel for the star formation, even if indi- 
rectly. Due to this fact, we now focus on the possible perturba- 
tions to the star formation activity that might be induced by the 
removal of the neutral hydrogen in the environment of a rich 
cluster of galaxies. The most suitable parameter to quantify the 
H I removal in a galaxy is the H I deficiency, computed according 
to lHavnes & Giovaneliild 19841) as the logarithmic difference be- 
tween Mm of a reference sample of isolated galaxies and Mm 
actually observed in individual objects: defm - log Mm ref - 
log Mm obs- log Mm ref nas been found to be linearly correlated 
to the galaxy linear diameter d as log Mm ref - a+b log d, where 
d (kpc) is determined at the 25''' B-band isophote, and a and b 
are weak functions of t he Hubble typ e . Here we adopt the up- 
dated coefficients from Sola nes et alj (Il996l) . In the following 
analysis we choose to consider a galaxy HI deficient if the ratio 



between the measured and the expected HI mass is > 2.5, i.e. if 
defm ^ 0.4. Therefore, our sample is composed of 54% non de- 
ficient galaxies that we consider unperturbed and 46% deficient 
galaxies that suffer from significant HI removal. 



3.2.1. Global abundances 

Firstly we focus on the relationship between the global SFR, H2 
mass and stellar mass as a function of the deficiency. Therefore, 
for each galaxy we compute a mean color along the radial profile 
between three pairs of quantities: sfr-sta, sfr-Hi and H^-sta, 
averaged over the entire disk. Since the color profiles typically 
exhibit a radial gradient that in some cases can be significant, 
considering a mean value for each galaxy may sometimes be 
a too crude assumption. We quantify this effect calculating the 
standard deviations <x and we obtain: 

0.07 < cr < 0.40 with a mean of cr = 0.20 dex for sfr - sta, 
0.05 < cr < 0.40 with a mean of & = 0.25 dex for sfr - H2 and 
0.03 < cr < 0.60 with a mean of cr — 0.31 dex for H2 - sta; 
these deviations, plotted in FigJU provide an estimate of the 
uncertainty in addition to the errors in the profile calibrations. 
To clarify the physical meaning of these colors we notice that 
sfr - sta is similar to the Ha equivalent width. In fact the sfr - 
sta color and the Ha equivalen t width (E.W.) (retriev e d from 
GOLDMine or co ll ected from [K oopman n & Kennevl (|2004); 
iRomanishin] (1 19901) : iKennicutt & Kent! d 1983b ) are found to be 
tightly correlated (r = 0.84). Similarly the H2 - sta color quanti- 
fies the H2 emission, normalized to the stellar continuum; finally 
sfr - H2 can be considered a tracer of the efficiency at which 
the molecular gas is converted into new stars. It is known that 
galaxies obey empirical scaling laws according to their dynam- 
ical mass and therefore to their H-band lumino sity, owing to 
the tig ht correlation between these two quantities dGavazzi et al.l 
1996a). As the luminosity increase s, the HI mass to luminosity 
ratio Mhi/Lh (Gavazzi et al. 2008), the H2 mass to luminosity 
ratio M ft /Lfl C Boselli et al. 1997) and the Ha equivalent width 
dGavazzi et al.l Il996al) decrease: in other words, the smallest 
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Fig. 4. The relationship between defm and color residuals (7/ 2 - sta)g (right panel) and (sfr - sta)n (left panel). Correlation 
coefficients are respectively r = 0.60 and r = 0.65. Uncertainties are discussed in the text. 



galaxies are, relative to their mass, the richest in gas and the most 
active in star formation. Before studying if the above three sets of 
colors display a residual dependence on defm independently of 
their mass, we need to correct for the primary luminosity depen- 
dence. Unperturbed galaxies in our sample (defm < 0.4) obey 
the following relationships, derived using a least square fit: 

H 2 - sta = -0.25 x log Lh + 1.55 (rms = 0.30dex) (10) 

sfr - sta = -0.26 x log L H + 1 .95 (rms = 0.32dex) . (11) 

In Figure|4]we plot, as a function of the HI deficiency, the residu- 
als (H2-sta)ft (left panel) and (sfr-sta)^ (right panel) computed 
using the previous equations. Inspecting the left panel, we can 
see a significant correlation (r = 0.65) between the (sfr - sta)n 
and the deficiency; since environmental processes do not sig- 
nificantly affect the stellar distribution in massive galaxies, the 
meaning of this correlation is that HI removal causes an increas- 
ing reduction in the SFR with increasing deficiency. In the right 
panel a weaker but still significant correlation (r = 0.60) exists 
between (H2 - sta)u and the defm, showing that molecular hy- 
drogen, even if well bounded in the galaxy potential well, can be 
perturbed by the HI removal, at least indirectly. Conversely we 
find (not shown) that sfr - H2 and defm are uncorrected quan- 
tities (r = 0.20), suggesting that the efficiency of molecular gas 
to star conversion does not depend significantly on the amount 
of HI removal. 

3.2.2. Radial behavior 

After having shown that both the SFR and the molecular abun- 
dance are reduced in deficient galaxies on a global scale, a new 
question arises: is this reduction mainly due to a global quench- 
ing or it does take place primarily in the external disk, producing 
a truncation in the star forming disk? This issue has been ad- 
dressed, among others, in a series of studies by Koopmann et al. 
They find that Virgo galaxies have their Ha disk truncated, while 
the case of starvation that produces a glo bal reduction in star 
formation across the entire disk is rare dKoopmann & Kennevl 



l2004h . To assess this result in our sample and to verify if similar 
conclusions hold also for the molecular gas, we study the radial 
behavior of the color profiles. We normalize the radius of each 
galaxy using r 2 5 from Tableland we compute the mean colors in 
bins of r/r 2 5, averaging over all galaxies in each bin. In Figure|5] 
we show a comparison between the mean sfr- sta profile for de- 
ficient (dashed lines) and non deficient (continuous lines) galax- 
ies; from left to right, the deficiency threshold increases, i.e. in 
the left panel deficient galaxies are considered with defm > 0.4, 
in the middle those with defm t. 0.65 and in the right panel 
the ones with defm > 0.9. Conversely, non deficient galaxies 
have defm < 0.4 in all panels. The first panel shows a gap of 
0.4/0.5 dex between the sfr - sta color for non deficient and 
deficient galaxies, suggesting that the SFR in deficient galaxies 
is on average reduced. Increasing the deficiency threshold (mid- 
dle panel) the gap increases up to 0.6 dex, while for the highest 
deficient galaxies defm ^ 0.9 (right panel) the gap increases up 
to 0.8 dex. Moreover, there are no galaxies with Ha emission 
beyond r/r 2 5 = 0.8; therefore, for extremely deficient galaxies, 
besides the global quenching of the star formation rate, there 
is an evident truncation effect. This result suggest that the cor- 
relation found between the global sfr - sta and the defm is 
mainly driven by quenching, while truncation is important only 
for highly deficient galaxies. For i/ 2 - sta, an effect similar to 
the one described for the sfr - sta is present in Figure[6] In fact, 
in the left panel the gap between non deficient and defm > 0.4 
galaxies is 0.4/0.6 dex, it increases to 0.8 dex when defm ^ 0.65 
(middle panel); for defm ^ 0.9 (right panel) the gap increases 
significantly with increasing radii and there is no emission be- 
yond r/r 2 5 = 1. Once again, this finding is consistent with the 
fact that the correlation found between the global i/ 2 - sta and 
the defm is mainly driven by quenching, while truncation is im- 
portant for highly deficient galaxies only. Conversely we find no 
significant difference between non deficient and deficient galax- 
ies as far as the radial variation of the sfr - H2 efficiency is 
concerned, being on average within + 0.2 dex. The existence of 
a significant correlation between the extent o f the HI disk and the 
de fm parameter is a well known result (e.g. lCavatte et all l994) 
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Fig. 5. Comparison between the mean sfr - sta profile for deficient (dashed lines) and non deficient (continuous lines) galaxies; 
from left to right, deficiency threshold is set to be defm > 0.4, defm > 0.65 and defm > 0.9. Non deficient galaxies are always 
considered for de f HI < 0.4. 
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Fig. 6. Comparison between the mean H2 - sta profile for deficient (dashed lines) and non deficient (continuous lines) galaxies; 
from left to right, deficiency threshold is set to be defm > 0.4, defm > 0.65 and defm > 0.9. Non deficient galaxies are always 
considered for defm < 0.4. 



implying that environmental effects, and in particular ram pres- 
sure stripping, perturb the ISM distribution outside-in, as also 
found in the simulations (lAbadi et al.l 11999b . In this work we 
consider the HI isophotal radius, taken arbitrarily at 5 x 10 19 
cm -2 , as a suitable measure of the HI disk. This threshold well 
matches the sensitivity of most of our HI data, except for NGC 
3521, 4192, 5055, 5236 and 5457 where the HI column density 
stays above the threshold at all radii, but their exponential decay 
can be confidently extrapolated outside; only for NGC 4402 and 
NGC 4579, whose HI data are limited to the very inner regions, 
would any extrapolation be meaningless and we exclude them. 
In Figure [7] we show the tight correlation (r = 0.81) between 
the isophotal HI radii normalized to T25 and defm, allowing us 
to quantify that at very large deficiencies (defm ^ 0.9) it occurs 
that the extent of the HI disk becomes smaller than the optical 
radius. Is there something fundamental in the transition between 
galaxies that have their HI disk smaller or bigger than the opti- 
cal radius? Does this transition also govern the difference in the 
SFR and H2 radial profiles between extremely deficient and nor- 
mal galaxies? To test these hypotheses, we further investigate if 
the outer parts of the SFR and H2 radial profiles exhibit a differ- 
ent behavior as a function of the deficiency. We derive for each 



galaxy the radial gradient by fitting the color profiles (sfr - sta 
and H2 - sta) in the outermost regions (r/r25 > 0.50). The fit is 
performed with a maximum likelihood method, assuming a con- 
stant uncertainty along the profiles and using an exponential, a 
reasonable model for the H2 - sta since both the molecular and 
the stellar component typically behave as exponential disks. The 
issue is more complex for the sfr - sta since the Ha emission 
tends to be rather erratic on small scales; however, after a 15 arc- 
sec smoothing even the Ha profiles can be fitted with an expo- 
nential function. The results of this analysis are summarized in 
Figure [8] where we plot against defm the mean radial gradients 
for the outer part of the sfr - sta (top panel) and H2 - sta (bot- 
tom panel) color profiles; the mean values are derived averag- 
ing in four bins of deficiency (0-0.4, 0.4-0.65, 0.65-0.90, >0.90) 
and the error bars refer to the uncertainties on these averages. 
The first panel of Figure [8] shows that the mean gradients of the 
sfr - sta color profiles are similar for non deficient and moder- 
ately deficient galaxies (up to defm — 0.90), while when consid- 
ering highly deficient galaxies the discrepancy becomes signifi- 
cant at lcr. Finding steeper sfr - sta profiles for highly deficient 
galaxies implies that when the defm > 0.9 and thus the HI disk 
is reduced inside the optical radius, the SFR sharply decreases in 
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Fig. 7. The relationship between defm and the HI isophotal ra- 
dius at 5 x 10 19 cm 2 , normalized to the i25- Correlation coeffi- 
cient is r — 0.81. The dashed lines emphasize that r HI >i25 for 
defm > 0.9. 

the outer part of the disk, consistently with the truncation effect. 
The second panel shows that the discrepancy between the mean 
gradients of the H2 - sta color profiles in non deficient or mod- 
erately deficient galaxies and extremely deficient galaxies is less 
sharp. In fact, due to the great dispersion in the 0.65-0.90 and 
>0.90 bins, the difference between highly deficient galaxies and 
moderate deficient galaxies is less marked; however, it is evident 
that the steepest H2 - sta are found in galaxies with defm > 0.9, 
once again supporting the idea that also the H2 profiles sharply 
decrease in the outer part of the highly deficient disks, suggest- 
ing that the truncation may also affect the H2 distribution. 

4. Discussion 

4.1. The Schmidt law 

The study of the local Schmidt law addresses one of the main 
open questions about the star formation activity in spiral galax- 
ies: is the process of star formation related to the atomic or the 
molecular component? Clearly, the formation of stars takes place 
inside the GMCs, in cold and dense molecular regions on a scale 
of few parsecs; however, while studying the molecular distribu- 
tion in our sample at a resolution of 15", our observations are 
sensitive to the diffuse component on typical scales of 1 kpc; 
this molecular medium can sometimes be not self gravitating 
and therefore it may not take part to the star formation activ- 
ity (lElme greed |2002|) . From our analysis of the local Schmidt 
law we find that there is a good correlation between the SFR 
and the H2, that eventually becomes poorer in the outer part 
of the disks (see Fig. [2] left). Since the bulk of the star for- 
mation activity takes place in the inner disks, it can be con- 
cluded that the Schmidt law for the molecular component is a 
good model of the SFR in galaxies. Therefore it seems plausible 
that the molecular medium on the kpc scale traces the self grav- 
itating molecular clouds which form at much smaller scales and 
sustain the star formation. However we found that the Schmidt 
law computed for the total gas, i.e. including the atomic hy- 
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Fig. 8. Mean radial gradients for the outer part of the sfr - sta 
(top panel) and H2 - sta (bottom panel) colour profiles, in four 
bins of deficiency (0-0.4, 0.4-0.65, 0.65-0.90, >0.90). Error bars 
are the uncertainties on the mean values in each bin. 



drogen, is significantly improved, particularly in the outer disks 
(see Fig. [2l right). Th is result supports the theoretical argument 
by lElmegreenl d2002l) that some clouds can be in a molecular 
state without being strongly self-gravitating and, on the con- 
trary, there might be instabilities in the HI component capable 
of triggering the star formation via direct collapse of HI; there- 
fore it is reasonable to accept that the star formation rate scales 
with the total column density of gas in all forms, including a pos- 
sible contribution from photodissociat ion of H2 molecule s. We 
thus disagree with the conclusion of Wong & Blitz (2002) who 
state that the correlation in the Schmidt law is driven entirely by 
the molecular component. The fact that we find a mean index 
n — 1 .4 for the local Schmidt law, similar to the one originally 
proposed by Kennicutt using integrated quantities, seems to im- 
ply that this relation is scale invariant, i.e. it is able to describe 
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the process of star formation on scales ranging from 1 kpc to the 
entire galaxy disk. 

4.2. Environmental effects 

The HI is the primary fuel for the star formation in galaxies, 
both directly and indirectly. We have shown that at the disk 
periphery the HI might directly contribute to the star forma- 
tion by collapsing into molecular clouds and indirectly by pro- 
viding the constituent of the diffuse H2 medium, composed of 
molec ules which form on dust gra ins inside the stellar potential 
wells (Blitz & Rosolowsky 2004). Therefore, it is not surpris- 
ing that the HI removal caused by the environmental processes 
may affect the star formation activity in cluster galaxies; we have 
shown that HI deficient galaxies have their SFR and H2 abun- 
dances reduced with respect to normal galaxies. Given the fact 
that we find no significant dependence of the SF efficiency on 
the environmental conditions, the observed trend of decreasing 
SFR and H2 density in deficient galaxie s is consistent with mod- 
els of galaxy starvation, as propo sed bv lLarson et al.l d!980t) and 
elaborated bv lBekki et al.l (12002b . according to which the effect 
of reducing the HI halo halts the replenishment of H2, leading 
to a quenching of the star formation rate. This model implies a 
continuous HI inflow from outside, as observed in some nearby 
galaxies, although at a rate th at seems insufficien t to sustain the 
current rate of star formation (ISancisi et al.ll2008l) . When galax- 
ies become extremely HI deficient, i.e. when the hostile envi- 
ronment has produced the ablation of their neutral gas up to a 
factor of ten and the HI disk reduces its size within the optical 
radius, superimposed on starvation, significant truncation of the 
star formation appears to exist, i.e . the SF disk is sharply erode d 
outside-in. The disagreement with Koopma nn & Kennevl (l2004). 
who maintain that the reduction in the global star formation in 
the Virgo cluster is caused primarily by the spatial truncation of 
the star-forming disks, and not due to starvation, is only sug- 
gested. Their sample is dominated by extremely deficient spirals 
(among their deficient galaxies, 65% have defyi > 0.9 and 85% 
have defm > 0.8). Beside the agreement with previous stud- 
ies on the SFR in clusters, our analysis points out for the first 
time that the H2 content might be affected by environmental pro- 
cesses, even if indirectly as a consequence of HI depletion. 

5. Conclusions 

With the aim of exploring the relations between the process of 
star formation and the physical properties of the ISM in various 
galaxy environments, we collected state-of-the-art imaging ma- 
terial and maps for 28 massive spiral galaxies belonging to the 
Virgo cluster and to the local field. The observational material 
includes: images of the stellar continuum (taken in the red/NIR 
bands); Ha images of the young (< 4x 10 6 yrs) stars; radio maps 
at 2mm from the recent Nobeyama CO atlas of nearby spiral 
galaxies that combines good sensitivity to large-scale CO emis- 
sion with a 15 arcsec spatial resolution, and sensitive radio maps 
at 21 cm from the ongoing, yet unpublished, VIVA and THINGS 
surveys carried out at the VLA. Physical parameters have been 
derived homogeneously and carefully for the individual galaxies 
applying corrections based on measured quantities, rather than 
average values. We hope that the effort we put in the analysis 
of the data reduces the statistical limitations that arise from the 
paucity of imaging material currently available. 
The present analysis indicates that the bulk of the star formation 
in spiral galaxies is supported by a diffuse molecular medium 
which forms through the conversion of the atomic hydrogen, due 



to the pressure exerted by the stellar potential. However, at the 
edge of the star forming disks, the HI plays a more important 
role than previously believed in directly sustaining the star for- 
mation activity. When environmental processes cause significant 
removal of the outer part of the HI disks, galaxies suffer from 
starvation; the replenishment of the atomic gas inside the optical 
disks is reduced, leading to a depletion of the molecular com- 
ponent that is consumed during the star formation activity, thus 
causing the quenching of the star formation activity itself. When 
the HI removal is so severe that the HI disk shrinks inside the 
optical radius, there is also a truncation of both the molecular 
and the star forming disk. 
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